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ABSTRACT 

In order to examine the "giant impact hypothesis" for the Moon formation, we run the first grid-based, 
high-resolution hydrodynamic simulations for an impact between proto-Earth and a proto-planet. The spatial 
resolution for the impact-generated disk is greatly improved from previous particle-based simulations. This 
allows us to explore fine structures of a circumterrestrial debris disk and its long-term evolution. We find 
that in order to form a debris disk from which a lunar-sized satellite can be accumulated, the impact must 
result in a disk of mostly liquid or solid debris, where pressure is not effective, well before the accumulation 
process starts. If the debris is dominated by vapor gas, strong spiral shocks are generated, and therefore the 
circumterrestrial disk cannot survive more than several days. This suggests that there could be an appropriate 
mass range for terrestrial planets to harbor a large moon as a result of giant impacts, since vaporization during 
an impact depends on the impact energy. 

Subject headings: planets and satellites: formation - method: numerical 



1. INTRODUCTION 

In the current standard scenario of planet formation, the 
final stage of assemblage of terrestrial plane ts is col lisions 
among proto-plane ts of abo ut Mars mass JKokubo & Idal 
ll998tlKortenkamp et al.ll2000t) . It is widely accepted that 
at this final assemblage stage, the Moon is formed from 
the circumterrestrial debris disk generated by an off-set 
impact of the proto-Earth with a Mars-sized proto-planet, 
which is known as the "Giant Impact (GI) hypothesis" 
JHartman & Davisll975llCameron & Wardll9 76). The initial 
phase of the GI scenario, a giant impact and the formation of 
the circumterrestrial debris disk, has been studied by a series 
of hydrodynamic simulations (|Benz 1 Slattery., & Cameronl 
fl986[lCameronl2000t lCanup & Asph augl200 It iCanutteOoUT 
N-body simulations of the accumulation of the Moon from 
the debris disk, which is the last phase of GI scenario, re- 
vealed that a single large moon is formed just outside the 
Roche limit, at a distance of about thre e to four times 
the Earth's radius, within se veral months ( Ida et alJ 119971 
Kokubo, Makin o. & Idal2 000). It is believed that this GI sce- 
nario expla ins a numb er of mysteries concerning the origin of 
the Moon (Cameron 2000): Why is the Moon so large com- 
pared to satellites of other planets? Why is the Moon deficient 
in iron and volatiles compared to the Earth? Why does the 
Earth-Moon system have large angular momentum? 

Almost all hydrodynamic simulations of GI in past decades 
used the Smoothed Particle Hydrodynamic (SP H) method , 
which is a particle-based L agrangian scheme l lLucvl 119771 
Gingold & Monaghanl ll977l) . In the SPH method, the nu- 
merical accuracy, in other words, the resolution, is deter- 
mined by the number of SPH particles. Th e latest simula- 
tions JCanup & Asphaudl2Wl1Cmiupll200l used 10 4 - 10 5 
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SPH particles, which is an improvement of 1-2 orders of mag- 
nitude compared to simulations in the previous tw o decades 
jBenz. Slattery & Cameronl! 19861 ICameronll2000l) . It is ap- 
parent, however, that the SPH method is not the best scheme 
to simulate an impact between two proto-planets. Firstly, SPH 
is not suitable to deal with the strong shocks and shear motion 
produced by the off-set impact. Moreover, since it is basically 
the Lagrangian scheme, the current SPH does not have resolu- 
tions fine enough for diffuse regions. This is critical problem 
particularly for GI, because the debris disk consists of only a 
few % of the total mass. As a result, even for a simulation with 
10 SPH particles, only a few 10 3 SPH particles are used to 
represent the debris disk, and thus the fine structure of the disk 
is not resolved accurately. In the SPH formalism, the spatial 
resolution is determined by the 'kernel size', which is variable 
with the local density. For a diffuse circumterrestrial disk, this 
kernel size can be as large as the ra dius of the disk itself (see 
Figure 1 of Canup & Asphaug 2001 ). This is the main reason 
that evolution of the debris is foll owed for very sho rt period 
(< 1 day) after the impact (Canup & Asphaus 2001). There- 
fore the post-impact evolution of the debris, which is a key for 
the GI hypothesis, is not well understood. 

Alternatively, one can use grid-based Eulerian methods to 
overcome the above problems. Melosh & Kipp ( 1989) have 
written the only preliminary study on the GI by the grid-based 
method to date. However, they neglected self-gravity of the 
planets and debris, and the evolution is followed only for less 
than 1 hour around the impact. Therefore, it cannot be com- 
pared to the recent SPH simulations. 

In this paper, we present the first three-dimensional hydro- 
dynamic simulations of the giant impact followed by forma- 
tion of a debris disk, taking into account the self-gravity, using 
a high-accuracy Eulerian-grid scheme. Our aim is to clarify 
the long-term (more than 100 hours after the impact) evolu- 
tion of the debris with the highest numerical accuracy used in 
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the previous simulations for GI. We pay attention especially 
to effects of pressure in the debris on formation of the large 
moon, rather than exploring the large parameter space for the 
mass ratio or orbits of the proto-planets. The previous SPH 
simulations with 10 4 - 10 5 particles would be fine for study- 
ing the behavior of GI for the initial several hours, and they 
suggest that the fraction of mass orbiting the proto-Earth just 
after the impact is mainly determined by the orbital parame- 
ters. On the other hand, the long-term behavior of the debris, 
and therefore the final mass of the moon depends on the evo- 
lution of the debris disk. It is expected that the pressure in 
the low-density debris should affect the hydrodynamic nature 
of the disk, e.g., by generating shocks. Pressure in the debris 
is determined by the thermodynamics and the phase-change 
of the material after the i mpact, which is not fully understood 
for GI (Stevenson 1987). For example, if a substantial frac- 
tion of the material is in a hot vapor form after the impact, the 
thermal pressure dominates dynamics of the disk, but if liq- 
uid or solid is a major component of the debris, its behavior 
could be very different. Since there is no practical numerical 
codes to treat directly two-phase flow and phase-change, we 
run hydrodynamic simulations assuming two extreme equa- 
tions of state, which approximately represent vapor or liquid- 
dominated materials. 

The paper is organized as follows. In §2, we briefly de- 
scribe our numerical method, equations of state, and initial 
conditions. We show the numerical results in §3. We give 
conclusions and discuss an implication on a necessary condi- 
tion of formation of large satellites for the Earth-type planets 
in §4. 

2. METHOD AND MODELS 

2.1. Numerical Methods 

The hydrodynamic scheme used here is a standard Eule- 
rian method, which has been widely applied to various astro- 
physical problems invo lving strong shocks in s t ructures with 
a high density contrast (Wada & Norman 2001; Wada 2001). 
We solve the following conservation equations for mass, mo- 
mentum, and energy, and the Poisson equation numerically in 
three dimensions: 

dp 



dt 



+ V-(/w) = 0, 



(1) 



dv Vp 
— + (v ■ V)v + — + V$ se = 0, 
dt p 



djpE) 
dt 



+ W-[(pE + P )v] = Q, 

V 2 $ sg = 47rG/9, 



(2) 

(3) 
(4) 

where, p,p,v, and E are the density, pressure, velocity of the 
matter, and the specific total energy. The gravitational poten- 
tial of the matter is denoted by $ sg . 

We use AUSM (Ad vection Upstream Splitting Method) 
( iLiou & Ste phen 1993| to solve the hydrodynamic equations 
and MUSCL (Monotone Upstream-centered Schemes for 
Conservation Laws) to achieve third-order spatial accuracy. 
Gravitational potential of the fluid on grid points is calcu- 
lated by solving eq. (@) using t he Fast Fourier Transfo rm with 
a convolution method (Hocknev & Eastwoo These 
are the standard techniques for hydrodynamics taking into ac- 
count self-gravity in grid-based schemes. More detailed de- 
scriptions on the numerical scheme and test calculations to 
guarantee the numerical accuracy are described in Wada & 
Norman (2001). 



We use an equally spaced Cartesian grid for the computa- 
tional box; 20x20x5 (model A and model B) in units of te, 
where is the radius of the proto-Earth. In order to see how 
the results are affected by a size of computational domain, we 
also run models in a larger domain, i.e. 40x40x 10(modelA' 
and model B'). The total number of grid points in the compu- 
tational box is 512 x 512 x 124 ~ 3.3 x 10 7 in all the models. 
The circumterrestrial region (i.e. r > r^) is represented by 
~ 10 7 grid points. 

2.2. Equation of State 

It is still an open question that what kind of equation of 
state (EOS) should be used in order to simulate complicated 
thermal and dynamical processes during GI and post-impact 
evolution. Conventionally two types of empirical EOS have 
been used with the SPH method: the T i llotson EOS (e.g. 
Canup & Asphaug 2001; Tillotsonl 1 19621 [Benz & Asphaug 
1 19991) and A NEOS and its extension (e.g. iMeloshl 120001; 
Canup 2004). ANEOS approximately describes mixed phase 
states (e.g., vapor and melt) by treating the different phases 
as separate components in temperature and pressure equilib- 
rium. Note that there are no numerical simulations for GI that 
can handle two-phase flow consistently, and that spatial reso- 
lution must be much smaller than the system size to describe 
the mixed phase states. As mentioned above, it is apparent 
that the spatial resolutions in previous SPH simulations are 
insufficient to describe mixed phase states, even if they use 
ANEOS (see also §3.3). 

One should also no t e that recen t SPH simulations 
(Canup & Asph aug||200Tl lCanur3 l2004) showed that the ini- 
tial behavior (i.e. a few dynamical times) of the impact does 
not strongly depend on choice of EOS for a given geometry of 
impact (i.e., mass ratio between proto-planets and the impact 
parameter). This suggests that the mass and the angular mo- 
mentum of the debris produced by the impact, which are the 
most important parameters for the accumulation process of a 
moon, are determined by a gravitational process, not by de- 
tailed thermal processes described by the conventional EOSs. 

We run a series of numerical experiments to follow a long- 
term (> 20 dynamical times) evolution of the debris disk. 
Instead of using the conventional EOSs, such as ANEOS 
or Tillotson EOS, we here assume two extreme EOSs: a 
polytrope-type EOS (EOS-1) and the same EOS, but with a 
zero-pressure cut-off (EOS-2). This is because assuming a 
simple EOS is more suitable to clarify essential physics be- 
hind the numerical results. EOS-1 is 



p = ( 7 - l)pE + C(n/3+l-j)p 



h/3+1 



(5) 



where E is the internal energy, and C is a constant. Through- 
out this paper, we assume the ratio of specific heats, 7 = 1 .01, 
and the constants C = 1, and n = 12. With EOS-1, the sys- 
tem behaves like an ideal gas or Polytropic gas at high or low 
temperature limits, respectively. EOS-2 is the same as EOS-1, 
but we introduce cut-off density (p c = 2.7 g cm" 3 ) and internal 
energy (E c = 4.72 x 10 10 erg g" 1 ), below which the pressure 
is forced to be zero in equations (j2ji and Q. Note that in 
the present problem, density in most of the debris is below 
Pq. These numbers are the sam e in the SPH sim ulations in 
which the Tillotson EOS is used ( iBenz & Asphau g8l999t) . We 
do not adopt the 'negati ve pressure' re gime assumed in the 
Tillotson EOS (Benz & Asphaug 1999), because it is appar- 
ently unphysical on the scale of GI. With this negative pres- 
sure regime, condensation in diffuse gas on any scales, e.g. 



Formation of Moon 



3 



the scale of the Moon, is allowed, even if self-gravity of the 
material does not work. 

The essential behavior of EOS-2 is similar to that of the 
Tillotson EOS. The cut-off for pressure phenomenologically 
represents a transition from vapor to liquid or solid particles. 
On the other hand, EOS-1 represents hot gaseous debris. The 
thermodynamical process during GI is complicated, for ex- 
ample, it is still unclear what fraction of the impactor is va- 
porized. This can be solved by a numerical scheme which can 
handle a two-phase (vapor and liquid) flow with elementary 
physics, but there is no practical code for the GI problem at 
this moment. We therefore perform numerical experiments 
for GI using hydrodynamic simulations with these two ex- 
treme EOSs to clarify the effect of pressure on post-impact 
evolution. 

2.3. Initial Conditions 

We follow lCanup & Asohaud J200H) and lCanudd200l for 
the orbital parameters of the impactor for which the most mas- 
sive satellite is expected. The masses of the proto-Earth and 
the impactor are assumed to be l.QM$ and 0.2M ffi , where M ffi 
is the Earth mass. The radii of the proto-Earth and proto- 
planet are te = 1.0 and 0.64r£, respectively. Note that no sig- 
nificant differences in the results for smaller impactors (e.g. 
O.lMg) were found in our simulations. The initial orbits of 
the impactor are assumed to be parabolic, and the angular mo- 
mentum is 0.86 Lg ra z, where L„„ is the angula r momentum 
for a grazing collision (Canup & Asphaug 2001 ). Initially the 
impactor is located at 4. Ore from the proto-Earth. 

3. RESULTS 

3.1. Disk Evolution and the Predicted Lunar Mass 

Figure 1 shows a typical time evolution of the giant impact 
with EOS-1 (m odel A). This model corr esponds to the 'late' 
impact model in Canup & Asphaug (2001). After the first im- 
pact (t ~ 1 hour), the disrupted impactor is reaccumulated to 
form a clump at t ~ 3 hours, which finally collides with the 
proto-Earth at t ~ 6 hours. During the second impact, the im- 
pactor is destroyed, and a dense part of the remnant spirals 
onto the proto-Earth (f ~ 10 hours), and a circumterrestrial 
debris disk is formed around t ~ 18 hours. It should be noted 
that many strong spiral shocks are generated in this process 
as seen in the density map (Figure 2) and azimuthal density 
profile (Figure 3). 

Figure 4 is the same model as model A (Figure 1), but with 
EOS-2 (model B). The initial behavior is similar to that of 
model A. The impactor is completely destroyed after the sec- 
ond impact as is the case in model A. However, the remnant 
is not diffuse, but more condensed (Fig. 4) and geometrically 
thin (see the edge-on views). This difference is reasonable, 
because with EOS-2, the internal pressure does not work in 
the low density gas. 

Next, we predict the satellite mass accumulated from the 
circumterrestrial disk in our models. Equation is an em- 
pirical formula for the satellite mass (M s ) as a function of 
the specific angular momentum of the circumterrestrial disk 
( /disk) derived from jV-body e xperiments of lunar accretion 
(Kokubo. Makino. & Ida 2003): 

-^-~1.9-^=-l.l-1.9^., (6) 

M disk VGM e a R M disk 

where M d j sk is the disk mass, a R is the Roche limit radius, and 
M esc is escaped mass, which is the mass lost during the ac- 
cumulation process. Using eq. (|6) and time-dependent mass 



and angular momentum taken from our numerical results as 
initial conditions for the accumulation process, we plot the 
time evolution of the predicted luna r mass in Figure 5. Here 
Mesc/Mftsk is assum ed to be 0.05 ( Kokub o". Makino. & Idal 
2000; Canup 2004). As seen in Fig. 5, the predicted lu- 
nar mass reaches a maximum of ~ 1 AMl for model A and 
~ \.2Mi for model B just after the second impact (t ~ 10 
hours). However, it should be noted that the predicted mass 
for the model A decreases very rapidly after the lunar mass 
reaches the maximum on a time-scale of 10 hours. The pre- 
dicted lunar mass becomes smaller than the current lunar mass 
Mi after t ~ 15 hours, and almost monotonically declines to 
0.3Ml at t ~ 65 hours. On the other hand, if we use EOS- 
2 which is relevant for a mostly liquid or solid material, the 
lunar mass increases again after t ~ 40 hours (model B). It 
reaches to ~ 0.6Mj,, and stays nearly constant until t ~ 90 
hours. 

In a fixed-grid Euler method, a choice of the spatial reso- 
lution and size of the computational box is a trade-off. If the 
computational volume is not large enough, the material that 
escapes from the computational domain in the initial phase of 
the impact cannot be followed, and this could reduce the fi- 
nal mass of the circumterrestrial disk. In order to check this 
effect, we also run model B', in which the computational vol- 
ume (40 x 40 x 10r|) is eight-times larger than that in model 
B (i.e. the grid size is twice larger than that in model B). The 
behavior is qualitatively similar with the result in model B, but 
the maximum of the predicted mass is about 1 .5 time larger 
(2Mi) than that in model B. This is because that a part of rem- 
nant 'arc' of the impactor is lost from the computational box 
in model B in the initial phase of the impact (see the snap- 
shots at t = 14.2 hours in Fig. 3). As a result, the predicted 
lunar mass in a quasi-steady state (t > 60 hours) is also 50% 
larger than that in model B. However, model B' also shows 
the lunar mass does not change after t ~ 50 hours, supporting 
the conclusion that the long-term evolution of the debris disk 
is followed correctly by the numerical resolution of model B' . 
Since the computational volume in model B' is large enough 
to cover more than 90% mass of the remnant arc formed by 
the impact, we expect that the predicted lunar mass cannot be 
much larger than Q.9Ml, even if we use a much larger compu- 
tational volume with the same spatial resolution in model B. 
Similarly, model A' covers 8 times larger volume than model 
A, and the peak lunar mass is 1.6 times larger than the one in 
model A. However, it also shows rapid decrease of the lunar 
mass, which is in contrast with the behavior in model B and 
B'. 

3.2. Angular Momentum Transfer by Spiral Shocks 

The reason why the predicted lunar mass cannot stay nearly 
constant after the impact in model A is clear. As seen in Fig- 
ures 2 and 3, spiral shocks are generated in the debris disk. 
Figure 3 shows that at least two clear shocks are present for 
the same radii (i.e. cf> = 0.2ir and 1.95 ir for r = 2.9rg and 
4> = 0.57T and 1.57T for r = 5. Org). The density jump at each 
shock is a factor of 4 to 16, suggesting that the effective Mach 
number in the model is about 2 to 4. In fact, a significant frac- 
tion of the debris is supersonic in model A. Figure 6, which 
is a frequency distribution of Mach number A4 in the debris 
at two different epochs, clearly shows that A4 > 6 (mass- 
weighted) or M. > 3 (volume-weighted) in majority of the 
debris at t = 12.3 hours. Even in a later stage, i.e. t = 57.1 
hours, the disk is still supersonic A4 >2. 

Due to these shocks in the gaseous disk, the inner massive 
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disk inside the Roche limit effectively falls onto the proto- 
Earth in a few rotational periods (see discussion below). This 
initial fall of the disk cannot be avoided in the GI process, if 
majority of the debris is in a hot gas phase. Since the colli- 
sion should be off-set, the generation of the spiral shocks is 
inevitable in a vapor gas disk with a high Mach number. 

The angular momentum transfer by the spiral shocks is 
more effective for stronger shocks with larger pitch angles 
(i.e., more open-spirals). We simply estimate it as follows. 
For isothermal shocks, the post-shock velocity is ~ vq/M 2 , 
where Vo is pre-shock velocity of the gas, perpendicular to the 
shock front. Therefore the stream line is bent by passing an 
obliqu e shock (see Fig. 12 and Appendix in Wada & Koda 
2004). In a rotating medium, this effect removes the angular 
momentum of a fluid element when it passes a standing spiral 
shock. The angular momentum change due to one passage of 
a spiral shock can be approximately estimated by 



. / sin 2 / , 

7/70 ra [iw +C0Sl 



1/2 



; - arctan 



.M 2 cos/ 



,(7) 



where i is the pitch angle of the spiral shock, jo and j 
are the angular mo mentum of pre- and post-shock flow 
( Wada lfe Kodall2"o"04l) . The sound velocity of SiCh gas with 
2000 K is 0.6 km s _1 for 100% vapor, and ^0.1 kms" 1 if th e 
mass ratio between vapor and liquid is 0.1 JStevensonll987l) . 
For strong shocks (A4 ^ 1), the angular momentum change, 
equation |7}, is j / jo ~ cos 2 i. Therefore, the angular momen- 
tum loss is ~ 12% per spiral shock with i = 20°. For weak 
shocks (e.g., M. = 2), equation @ suggests that the angular 
momentum loss per shock is about 5.6% for i = 20° and 1.4% 
for i = 10°. As shown in Fig. 2 and Fig. 3, several spiral 
shocks are generated in the debris disk, then a fluid element 
would lose ~ 10% of its angular momentum for one rotational 
period. This is large enough for the debris gas to lose most of 
their angular momentum in 100 hours after the impact. There- 
fore, the presence of spiral shocks in the debris disk is crucial 
for evolution of the disk, and as a result for the mass of satel- 
lites. 

3.3. Comparison with Previous SPH Results 

Remarkably, although we use methods, resolutions, and 
equations of state different from previous SPH simulations, 
the predicted mass just after the impact is consistent with 
recent SPH results: 1 .7 — 1.8Mf for the 'late impact' model 
(Canup & Asphaug 2001: ICanupll2004l) . This again implies 
that the early phase of the impact process (t < 10 hours) is 
dominated by gravity, but not by hydrodynamical and thermal 
processes. However, we find that behavior of the remnant and 
the predicted lunar mass in the late phase depend on EOS, 
especially whether pressure works in the disk. This was not 
pointed out by the previous SPH simulations. One obvious 
reason is that the spatial resolutions in the SPH results are 
not fine enough to follow long-term evolution of the debris 
until the effect of pressure bec omes evident. In the previous 
SPH simulations (Canup 2004), only ~ 10 3 particles are used 
to represent the debris. Suppose 1000 SPH particles are uni- 
formly distributed in a thin disk inside the Roche limit (3r^), 
the average separation between particles is ~ 0.2r^. In SPH, 
variable kernels are often used, and the kernel size is roughly 
the same as the spatial resolution. Usually an SPH particle 
is assumed to interact with 30-50 neighbor particles. There- 
fore, the kernel size in the previous SPH simulations, may 
be ~ in the disk, which is about 10-20 times larger than 



our grid size. Moreover, the debris has in fact an extended 
distribution, rather than a thin disk, and the SPH kernel size 
is much larger than the numbers estimated above. It is not 
straightforward to compare the spatial resolution between the 
two different schemes, but apparently that the spatial resolu- 
tion in our simulation is greatly improved from the previous 
ones. For example, in Figure 3, the azimuthal density pro- 
file is obtained by about 1000-2000 grid points. However, in 
the previous SPH simulations, only 10- 100 particles would be 
available for drawing this kind of density profile. 

The other reason why the previous studies did not notice the 
effect of shocks is that the conventional EOSs (ANEOS and 
Tillotson) are close to our EOS-2, in a sense that pressure is 
not effective in the diffuse debris. 

4. CONCLUSIONS AND DISCUSSION 

We run for the first time three-dimensional grid-base hy- 
drodynamic simulations of the giant impact between proto- 
Earth and a proto-planet taking into account self-gravity. We 
assume two types of the equation of state that phenomeno- 
logically represent hot gaseous material or liquid/solid mate- 
rial. Our numerical experiments suggest that in order to form 
a lunar-sized satellite from the circumterrestrial debris disk 
produced by the giant impact, the most fraction of the debris 
should not be in a pressure-dominated phase (e.g., hot vapor 
gas). Otherwise the subsequent disk evolution results in form- 
ing only a small satellite because of the fast angular momen- 
tum transfer associated with spiral shocks prior to the satel- 
lite accretion stage. The time-scale of the angular momentum 
transfer is an order of 10 days, therefore the accumulation to 
form a large satellite should be much faster than this. Yet our 
results may not rule out a possibility forming a satellite as 
large as the Moon from a pressure-dominated disk, if the disk 
mass is much larger than the current lunar mass. This might 
be the case for collisions with a massive impactor, but this 
causes another problem on a fraction of vapor in the debris 
(see discussion below). 

In agreement with the previous S PH simulations (e.g., 
Canup & Asphauel l2"o"Oll ICanupll20 04), we find that the pre- 
dicted lunar mass at t ~ 10 hours after the collision does not 
strongly depend on choice of EOS and numerical methods for 
the same orbital parameters. This means that the early phase 
of GI is dominated by a gravitational process, not by thermo- 
dynamical processes. On the other hand, the late phase of GI, 
i.e., evolution of the debris disk, is sensitive to EOS, espe- 
cially the pressure in the disk. In order to clarify this differ- 
ence, the spatial resolutions of the previous SPH simulations 
were apparently not fine enough. 

Our results give an important implication on a necessary 
condition for formation of large satellites for the Earth-type 
planets. The present results suggest that one of the important 
keys for GI scenario is the fraction of the vaporized debris af- 
ter the giant impact. If the kinetic energy of the collision is 
much larger than the latent heat of the major component of 
the proto-planet, most of the proto-planet could be vaporized 
by the impact. If this is the case, the disk evolution would 
not lead to formation of a large moon as explained above. On 
the other hand, if the impact velocity is slower than a criti- 
cal value, a large fraction of the material can be in a liquid or 
solid phase, and shocks do not dominate the angular momen- 
tum transfer in the debris disk. In this case, the post-impact 
evolution of the proto-planet could be similar to that in model 
B in our experiments, and therefore a single large moon could 
be formed. By assuming the impact is head-on and compar- 
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ing the latent heat of SiCh JStevensonl ll987) and the kinetic 
energy of the impact, we can roughly estimate the critical im- 
pact velocity as ~15 km s" 1 , which is slightly larger than the 
surface escape velocity of the Earth. In the grazing impact 
between proto-Earth and a Mars-sized proto-planet, shocks 
propagating in the planets involve vaporization, and its veloc- 
ity is probably comparable (or smaller) to this critical value, 
and therefore it would be natural to postulate that most of the 
debris mass is in a liquid phase. This argument leads to an 
interesting suggestion: if the impact velocity is larger than the 
critical velocity, in other words, the mass of the planet is larger 
than, say, a few Earth mass, the giant impact never results in 
forming a large satellite. 

Finally, one should recall that a correct phase-change even 
for the major component, such as SiC>2, during GI followed 
by for mation of a debris disk is still unknown (e.g. Stevenson 
1987). For example, the post-impact expansion could pro- 



duce more vapor in the debris material, and this may change 
the hydrodynamic property of the circumterrestrial disk. The 
final fate, namely whether the planet has satellites, and how 
large they are, depends on the nature of the phase-change. 
The thermo-dynamical process during catastrophic impact be- 
tween the proto-planets is still too complicated to be explored 
using current numerical techniques including the SPH and 
grid-based methods. This ultimately requires a self-consistent 
numerical scheme to simulate two-phase flow taking into ac- 
count the realistic thermodynamical processes. 
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Fig. 2.— 



Snapshot of the density field of model A at r = 12.3 hours. Strong spiral shocks in the debris are resolved. 



Formation of Moon 
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FIG. 3. — The azimuthal density profiles at r = 2.9 ye (upper dots) and 5.0 (lower dots) on z = plane for the same density snapshot in Fig. 2. The unit of 
density is po = 12.6 g cirT 3 . 





FIG. 4. — Same as in Figure 1, but for model B, in which zero-pressure EOS (EOS-2) is assumed. Following previous simulations with the Tillotson EOS 
i Tillotson 1962), we adopt the cut-off density [p c = 2.7 g cm~ 3 ) and internal energy (E c = 4.72 X 10 10 erg g~'), below which the pressure is forced to be zero in 
EOS-1. The EOS-2 represents a state where most of the mass is in a state of liquid, not in vapor. 
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FIG. 6. — Frequency distribution of Mach number in the debris (p < O.lpo) of model A. Thick lines are at ( = 12.3 hours. Thin lines are at t = 57.1 hours. 
Solid and dotted lines are mass-weighted and volume-weighted histograms, respectively. 



